function rho = dpropa(nd,rho,omsc,gamsc,nleja,zl,zdvd,cset,coef)
% Calculate the action of the expansion of the function of the
% Liouvillian in terms of Newton polinomials over the density matrix.
% raux auxiliary matrix.

% First term (f[z1]1_op rho).
      
        raux = zdvd(1)*rho;

% Rest of terms.

      for k=2:nleja

% Calculating (L - z_(k-1)1_op) rho.
   rho = newton(nd,rho,omsc,gamsc,zl(k-1),cset,coef);

          raux = raux + zdvd(k)*rho;
 
      end
 


% raux -> rho
        
            rho=raux;

      end
